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I . INTRODUCTION^ 

The far-field spectrum of a supersonic jet has a distinctly 

different spectrum than that of a subsonic jet. The spectrum of a 

supersonic jet diminishes with increasing frequency at a 1/f rate 

(where f is frequency) whereas a subsonic jet spectrum has a slope 
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proportional to 1/f . It was proposed that the difference could be 
attributed to a nonlinear acoustic rather than aerodynamic effect. 
Research on high intensity (finite amplitude) acoustic waves has 
shown that nonlinear distortion effects generally result in a shift 
of energy to higher frequencies. The higher intensities associated 
with supersonic jets would therefore indicate that high frequency 
enhancement of the spectrum should occur, resulting in the differences 
observed between subsonic and supersonic jets. 

Research on high intensity acoustic phenomena has been largely 
restricted to initially sinusoidal waves. The wave distortion effect 
in this case for both plane and spherical waves has been extensively 
treated theoretically and experimentally. 

The plane wave theory for initially sinusoidal propagation has 
been sufficiently verfied by experiment and the Bessel -Fubini solution 
(References 1 and 2) suffices for a large number of plane wave finite 
amplitude situations. This solution predicts that an initially sinu- 
soidal wave develops into a saw-tooth waveform or acoustic shock wave, 
thus the shift in energy from the sinusoidal fundamental wave into 
harmonics of higher frequency. 


Substantial portions of sections I and II of this report were presented 
in our renewal proposal and at the SHOCK NOISE WORKSHOP at NASA Lewis 
Research Center on December 5, 1980. 



Initially sinusoidal spherically diverging high amplitude waves 
behave in a similar manner; however, higher intensities are required 
for spherical waves to form shocks since the waves are diminishing in 
amplitude due to spherical divergence. The theory for the initially 
sinusoidal spherical case (reference 3) has also been verified exper- 
imentally (references 4-5). In both cases, for sufficiently high 
amplitude, the wave distorts into a stable saw-tooth form with a 1/f 
spectrum which, resembles the spectrum observed for supersonic Jet 
noise. 

The propagation of two superimposed initially sinusoidal high 
intensity waves of two frequencies has also been extensively investi- 
gated. The resxilting production of a subharmonic at the difference 
between the two frequencies was predicted theoretically by Westervelt 
(reference 6) and verified experimentally for propagation in water 
(reference 7). This phenomena, referred to as the ’’Parametric Effect”, 
is analogous to intermodulation distortion in electronic circuits and 
produces sum as well as difference frequencies. 

The finite amplitude effects for a single and dual frequency 
initially sinusoidal waves can be represented in terms of Fourier 
series for the sound pressure, p, as follows: 
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For the dual frequency case, there is both high and low frequency 
enhancement of the spectrum. The amplitude factors of the various 
harmonics for both the single and dual frequency cases are functions 
of the initial amplitudes and frequencies, the distance from the source 
z, the sound velocity c and the medium nonlinearity parameter which 
for the case of a gas is the ratio of specific heats. 

The dual frequency case has been observed in air by Bennett and 
Blackstock (reference 8) and their results are in good agreement 



with Che cheorecical CreaCmenc of Fenlon (reference 9) . The power 
llmltatlona of the sources used by Bennett and Blackstock were 
such that the difference frequency component was at least 50 dB 
lower than the primary levels generated at and They did 

not give results for uhe levels at u) . 
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It is reasonable to assume that for the case of propagation 
of high intensity broadband random noise that both the progressive 
wave steepening effect observed for single frequencies and the 
intermodulatlon effect of dual frequency propagation are occurring 
resulting in both low and high frequency spectrum enhancement. 

The propagation of high intensity random noise has received 
considerably less attention that the propagation of Initially 
sinusoidal waves; however, several significant studies have been 
conducted. Pernet and Payne (reference 10) demonstrated that for 
a relatively narrow band of high intensity random noise, harmonic 
bands are generated somewhat analogous to the periodic case. This 
study was limited to plane waves. In another plane wave study, 
Pestorius and Blackstock (reference 11) found that for sufficiently 
high levels of broad band noise that shocks are formed and that 
there is energy conversion to both the high and low frequency 
directions; however, the predominant conversion effect is in the 
high frequencies with a distinct change in the slope of the spectrum. 

A study of the propagation of spherically diverging high 
amplitude random noise has been conducted by Webster and Blackstock 
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(reference 12) . Sound pressure levels of 121 to 145 dB at one meter 
from the source were generated using a conventional electroacoustic 
transmitter over a frequency range of 2-10 kHz. Although the acoustic 
levels were not sufficient to produce acoustic shocks, strong 
nonlinear behavior was nevertheless observed. Nonlinear high 
frequency effects were observed over the entire source to receiver 
propagation path and Webster and Blackstock conduce that the 
nonlinearly generated high frequency noise never becomes small 
signal in behavior regardless of distance. 

Webster and Blackstock have also investigated the interaction 
of broad band noise with a sinusoidal signal for the plane wave 
case (reference 13) and have shown that the broad band noise is 
drastically altered in both spectrum shape and level by the presence 
of the periodic wave. This effect has, however, not been investigated 
for the case of spherical propagation. If this strong Interaction 
does occur for the spherical case it would imply that the charac- 
teristics of noise radiation from the exhaust of a jet engine could 
be altered by the presence of periodic waves originating from the 
compressor, turbine or oscillations due to instabilities. 

In this current research program the work of Webster and 
Blackstock on free-field propagation of finite amplitude noise 
(reference 12) is being extended to considerably high source 
intensity levels more representative of an actual supersonic jet. 
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A 10,000 acoustic watt source has been installed in an anechoic 
chamber. The levels generated by the source are sufficient such 
that acoustic shocks are readily obsei"ved. In addition, dual frequency 
excitation of the source has produced a strong parametric effect with 
a difference frequency comparable in level to the primary frequencies. 
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II. STATUS OF EXPERIMENTAL PROGRAM 

A Ling Model EPT-200 high intensity acoustic source has been 
installed in the Noise Control Lauoratory anechoic chamber. The 
associated operating system consisting of a 3.5 kw power amplifier, 
field supply and cooling system, and spectrum shaping equipment has 
been installed in the Laboratory Control room. An existing compressed 
air control system has been modified to provide the necessary air flow 
and pressure requirements of the source. The source is electrically 
modulated and can therefore generate both broad band and periodic 
waveforms. An exponential horn is coupled to the source to provide 
efficient acoustic power transfer to the free-field environment of 
the anechoic chamber. Photographs of the transducer installation, 
operating system, air education system and analysis equipment are shown 
in figures 1 and 4. The acoustic field is measured with a 1/4" B 5 K 
microphone mounted on a linear traversing mechanism. The results to 
date have been obt.iined by manually positioning the microphone at 
various distances from the source. For future experiments a radial 
traversing mechanism will also be installed and computer controlled 
stepping motors will be used for both linear and radial positioning of 
the microphone. In addition, controls for the air supply system are 
currently being computerized. 

The data analysis system consists of a Nicolet Model 660A dual 
channel Fast Fourier Transform Analyzer with a dedicated dual disk 
storage unit. The results are plotted with Tektronix digital plotter. 
The analyzer is interfaced with a Commodore 32K microcomputer which 
also has dual disc storage. 
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Preliminary results have been obtained for source excitation 
with band limited noise, sine waves, and dual sine waves. Typical 
results tor band limited noise are shown in figures S through 11. 

The spectrum shaper was adjusted for equal output at the 630, 800, 

1000, and 1250 Hz one-third octave bands with all other bands off. 

The w.iveform of the electrical input to the power amplifier is shown 
in figure 5. The total time span of this record and of figures 7, 9, 

and 11 is .008 seconds. The spectrum of this signal is shown in 

figure 6. The vertical scales on all spectra shown is 10 dB per 

division. The acoustic output at 3-1/2 teet from the source is shown 

in figures 7 and 8. The sound pressure level at this position is 
158 dB. ■f'he wave steepening effect is quite evident and the enhance- 
ment of high frequencies in the spectrum is present. The acoustic 
output at 7 feet is shown in figures 9 and 10. The waveform is steepep 
than at 5-1/2 feet and the high frequency content continues to increase 
with a change in slope of the spectnum. Continued steepening of the 
waveform is illustrated in figure 11 for a 10 foot microphone position. 

Spectra for a single 1/3 octave band of noise are shown in 
figures 12 through IS. A logarithmic frequency scale is used in this 
case. The electrical input spectrum is shown in figure 12. The spectra 
at .5 feet and 7,5 feet are shown in figures 13 and 14 and the waveform 
at ".S feet is shown in figure 15. The spectrum shift and steep wave- 
form are well illustrated 

The waveform and spectrum for a 1 kHz sine wave input are shown 
in figures lb and 17 measured at 4 ieet with a sound pressure level 
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of 158 dB. 

Dual frequency excitation is* shown in figures 18 and 19 for 
input frequencies of 800 and 1000 Hi. In this case the difference 
frequency component generated is higher than the components at 800 


and 1000 Hz. 










21 


III. STATUS OF THEORETICAL INVESTIGATION 

This study is largely based on an application of Burger's 
equation developed by McKendree and Fenton (reference 13). Burger's 
equation is a second order form of the wave equation. It describes 
suitably the propagation of sound in a non-linear medium. 

Burger's equation may be written as 

i£ _ o 

3a Pst r 3^ 

where for spherical waves 

r P' 

Pmax 

g = 3 e k Tf, In r 
^0 

T = 0) (t - — ) 

‘io 

r - ^ 

‘ Otw 

and 

ka^ 

r,j : Rayleigh distance = for a circular piston with radius a 
3 = Y = ^ * 1.4 in air 

Cy 

c = peak Mach number 

0 ^ = viscous absorption/wavelength 

p' i peak pressure 

r = radial distance 

Cq = phase velocity 


V 



22 


The term k ■ ^ has presented a problem. In the development of 
solutions to Burger's equation, this term has been considered arbitrary. 
Indeed, it is arbitrary for the plane wave solution. However, in the 
spherical wave solution it can play a very important role. If a shock 
will form, it will form at a » 1. To minimize the distance r this 
corresponds to k must be defined by 


k* 


4 

Bea' 


for a circular piston source. Due to spherical spreading, however, a 
shock may never form. 

Returning now to equation (1) . Assioming that terms can be separated 
we may apply 


lE - 

9a 
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over a small increment Aa and then apply 


( 2 ) 
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(3) 


over the same internal Aa. 

Equation (2) is simply the inviscid case of Burger's equation and 
may be solved using a finite difference solution. It then becomes after 
slight manipulation 

p(a,T) * p(o,T) + op(a,T) (4) 

where a is now used in place of Aa, a small increment. 

Since op(a,T) is small, and p is a well behaved function in a the 
terms in equation (4) may be approximated by the first two terms ct a 
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Taylor's series. Therefore, 

p(o,T) • p(o,T) ♦ op(o,T) || (5) 

Equation (5) is the basis of the present computer model. Viscous 
losses have not yet been included but may be using equation (3) over 
each increment of a in equation (5) . 

The model was first tested using an initially sinusoidal input 
waveform as shown in figure 20. Figures 21-23 show the waveform 
steepening and shoci, formation. The next step in this theoretical 
study will be to use actual experimentally measured random initial wave 
forms and compare the waveform and spectrum alteration predicted by the 
theory with the experimental data. 
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